Boundary-induced violation of the Dirac fermion parity and its signatures 
in local and global tunneling spectra of graphene 
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Extended defects in graphene, such as hnear edges, break the translational invariance and can 
also have an impact on the symmetries specific to massless Dirac-like quasiparticles in this material. 
The paper examines the consequences of a broken Dirac fermion parity in the framework of the 
effective boundary conditions varying from the Berry-Mondragon mass confinement to a zigzag 
edge. The parity breaking reflects the structural sublattice asymmetry of zigzag-type edges and 
is closely related to the previously predicted time-reversal symmetric edge states. We calculate 
the local and global densities of the edge states and show that they carry a specific polarization, 
resembling, to some extent, that of spin-polarized materials. The lack of the parity leads to a 
nonanalytical particle-hole asymmetry in the edge-state properties. We use our findings to interpret 
recently observed tunneling spectra in zigzag-terminated graphene. We also propose a graphene- 
based tunneling device where the particle-hole asymmetric edge states result in a strongly nonlinear 
conductance- voltage characteristics, which could be used to manipulate the tunneling transport. 
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PACS numbers: 73.20.At,73.22.Gk,73.63.Bd 



I. INTRODUCTION 

In condensed matter systems with nodal fermionic 
spectra, quantum description of low-energy excitations 
can resemble that of the ultrarelativistic electron. The 
crystal space group imposes a generic restriction on such 
quasiparticles known as fermion doubling: they come in 
pairs of opposite chirality species that can be mapped to 
the conventional "right-handed" (RH) and " left-handed" 
(LH) fermions of the Dirac theoryii The most recently 
studied examples are graphene, where two distinct Fermi 
points in the Brillouin zone give rise to both chiral 
speciesj^i^ and 2D HgTe quantum wells where spin-orbit 
coupling effectively results in a pair of the RH and LH 
fermions at low energies i^i^ The fermion doubling brings 
the symmetry with respect to the exchange of the chiral- 
ities, RH ^ LH, related to the parity symmetry of the 
Dirac equation.^ It is of considerable interest to inves- 
tigate the consequences of the violation of such a sym- 
metry, since they could be observable in materials where 
quasiparticles imitate Dirac electrons. Besides, the Dirac 
fermion parity is distinct from other discrete symmetries 
(e.g. time- reversal invariance), and, therefore, through 
its violation one could achieve additional control over 
electronic properties of the material. 

It has been noticed^i^ that discrete symmetries of a 
Dirac fermion system can be broken, along with the 
translational invariance, by the boundaries of the sys- 
tem. In this paper we focus on the parity violation due 
to such a boundary effect and suggest how to detect and, 
possibly, use it in electronic devices. 

Our main assumption is that the boundary does not 
cause scattering between the opposite chirality quasipar- 
ticles. To model this we use effective boundary con- 
ditions,—!^ interpolating between the infinite mass con- 
finemenlj^ and the zigzag graphene boundary.— The 
parity breaking occurs as long as the boundary devi- 



ates from the infinite mass confinement toward the zigzag 
edge. This is due to the structural sublattice asymme- 
try: a zigzag-type crystal face has unequal numbers of 
sites from the two sublattices of the honeycomb struc- 
ture.— More generally, the parity in this continuum 
model can cannot be preserved simultaneously with the 
time-reversal invariance near the edge of the system. 
This is closely related to the existence of time-reversal 
symmetric, propagating edge statesii which have become 
a topic of vigorous graphene-related research (see, e.g., 
Refs. 11114,15, 16, 17, 18, 19, 20, 9^21). The time-reversal 
symmetry requires that the edge states from the differ- 
ent valleys propagate in the opposite directions, forming 
a Kramers pair at a given energy. As a result, even in 
the absence of the intervalley scattering the problem does 
not reduce to a single valley, in the sense that the valley 
contributions to observables are not identical. For the 
low-energy states (imitating the RH and LH fermions), 
this implies broken Dirac fermion parity. 

We intend to demonstrate several new properties of the 
broken-parity edge states: 

(i) nonanalytic particle-hole asymmetry of local and 
global densities of states, 

(ii) time-reversal invariant pseudospin polarization 
(which in graphene is associated with the sublattice de- 
gree of freedom) , and 

(iii) asymmetric nonlinear bias-voltage dependence of 
the tunneling conductance. 

In view of the progress in experimental control over 
graphene edges;^ this material is particularly suitable to 
test our findings. Below we discuss in more detail the con- 
nection between our results and the ongoing graphene- 
related research. 

Our finding (i) can be tested by means of scanning 
tunnchng spectroscopy (STS) of the density of states 
(DOS). In fact, some STS experimentsA'^-^-^. have already 
reported a particle-hole asymmetric DOS with a peak 
at -20-50 meV for monoatomic zigzag graphene edges. 
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By contrast Li et ali^ have observed a symmetric lin- 
ear DOS in graphene bulk. In our model, the crossover 
from an asymmetric edge DOS to a symmetric bulk 
one follows naturally from the existence of the broken- 
parity Dirac fermion edge states. The results of experi- 
mentsi^ii^ could therefore provide some evidence for the 
Dirac fermion parity violation. Such an interpretation 
is also supported by the observation that the position 
of the DOS peak in the experiment of Niimi et a\J^ 
can be described very accurately by our model. This 
is achieved by taking into account not only the struc- 
tural asymmetry but also a potential-energy difference 
between the sublattices^ that generates weakly disper- 
sive edge states rather than the singular zero-energy 
band,^-'^ It was shown earlier-^ that the next-nearest- 
neighbor hopping on the honeycomb lattice could also 
result in a particle-hole asymmetric edge DOS. However, 
this symmetry-breaking mechanism results in the DOS 
peak at significantly larger energies, of the order of next- 
nearest-neighbor hopping energy w 300 meV. 

Another our result (ii) demonstrates that broken- 
parity edge states carry a time-reversal-invariant pseu- 
dospin polarization. This agrees with the general percep- 
tion that in graphene electronic properties and those aris- 
ing from the sublattice degree of freedom (pseudospin) 
are interrelated. However, there is a great deal of uncer- 
tainty as to how such a relation can be studied. The 
specific feature of our edge problem is that it is pos- 
sible to establish a one-to-one correspondence between 
the pseudospin polarization and the edge DOS. We sug- 
gest that the pseudospin polarization can be detected via 
measurements of the electric conductance in lateral tun- 
nel contacts between zigzag-terminated graphene and a 
suitably chosen metallic electrode. 

Our proposal is based on finding (iii) that the edge- 
state (i.e. polarization-dependent) contribution to the 
conductance is asymmetric with respect to the bias volt- 
age. Therefore, it can be separated from the symmetric 
contribution of the bulk graphene states. In addition, the 
edge-state tunnel conductance turns out to be strongly 
nonlinear: It exhibits kink-like switching as the sign of 
the voltage reverses. Such a behavior could serve as a 
prototype for the potentially useful electronic function- 
ality. 

The outline of the paper is as follows. In Sec. HIl we 
formulate the boundary problem for the Green's func- 
tion of the Dirac equation and discuss the role of the 
parity symmetry. In Sec. IIIII we analyze the local DOS 
and pseudospin polarization, and compare our results for 
the local DOS with the experimental data of Niimi et 
al.^^ Section [TVl describes the relation between the global 
edge DOS and pseudospin polarization and their tun- 
neling spectroscopy. The last section |V] summarizes our 
results and discusses their validity as well as possible ap- 
plications. 



II. DIRAC FERMIONS IN 2D SEMI-SPACE: 
BROKEN PARITY AND EDGE STATES 

Edge-state spectroscopy usually deals with isolated 
edges in large samples where finite-size effects are pre- 
sumably irrelevant fi^ii^iS^ We model this by considering 
a boundary problem for a Dirac fermion retarded Green's 
function in a 2D semi-space — oo<a;<cxD,0<j/<oo: 

(e/-i;75Sp)Ge(r,r') = 5(r-r'), (1) 

= (^Sn+ + ^Sn_) Ge|^^,, n^^ = 1, (2) 

with G|y^oo being finite. In Eq. ^ p — —ih{dx,dy,0), 
e and v are the 2D momentum operator, energy and ve- 
locity near a Fermi point. In view of the further analysis 
of the parity symmetry, equations ((!]) and ([2]) are both 
expressed in terms of the chirality, 7^ and effective spin, 
S through the Dirac matricesi^ 

75 = ^707172^3 ^ ^3 ^ ^0^ ^5^0^ = r° (g) cr,(3) 

7° = -r^ cr", 7 = ZT2®cr. (4) 

We introduce the two sets of Pauli matrices, cr^'^^s g^^^ 
T^'^'"^, and the corresponding unit matrices, a'^,T^ and 
/ = T° (g) tj". In graphene, cr^'^''^ represent the two sub- 
lattices of the honeycomb structure, while r^'^'-^ act in the 
valley space. The eigenstates and eigenvalues (r = ±1) of 
the hermitean matrix 7^ conventionally define the right- 
handed (RH, +) and left-handed (LH, — ) quasiparticles^. 
In Eqs. ^ and ([2]) they are described by the projected 
Green's functions i(/±7^)Ge. 

The boundary condition, Eq. ([2]) ensures vanishing of 
the particle current across the edge^ (i.e. no Klein tun- 
neling^"'). It is diagonal in chirality space with the RH 
and LH blocks parametrized by three-dimensional unit 
vectors, rir — n±, orthogonal to the boundary normal. 
In graphene, where t — ±1 label the valleys, the case of 
the zigzag edge corresponds tc^ 

n+ = -n_ = z, G, = 7's3Ge ly^g, (5) 

where z is the out-of-plane unit vector. This implies that 
one of the sublattice Green's functions must vanish at 
the edge, refiecting the structural sublattice asymmetry 
of the zigzag boundaryi^i The other nontrivial limit is 

n+=n_=x, Ge = ^'G,\y^^. (6) 

It is the infinite mass confinement of Berry and Mon- 
dragorti^ (x is the unit vector along the edge) . As shown 
m Ref. fgl, the intermediate case, when n± interpolate 
between JS]) and can be treated as a zigzag edge 
where the structural sublattice asymmetry coexists with 
a potential (energy) sublattice asymmetry that could re- 
sult from electron-electron interactions within atomic dis- 
tances near the edge.^^ This does not however exhaust 
the applicability of the boundary condition since it 
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can be derived from the only requirement that the Dirac 
particle current is zero in the normal direction at the 
edge.— 

Unless restricted to the infinite mass confinement case 
([6]), the boundary parameters are not identical, ^ n_, 
which is verified below on the basis of the time-reversal 
(T) symmetry [see, Eq. ()15p ]. Therefore, boundary con- 
dition ([2) (and, in particular, Eq. (O) explicitly contains 
the chirality, 7^. This violates the symmetry under the 
exchange of the RH and LH quasiparticles, 

7"G,7", (7) 

since it reverses the sign of 7^ (In contrast, S is even 
under such operation). On the other hand, the RH ^ 
LH exchange is involved in the parity transformation,— 

G,(r,r')^7"G,(-r,-r')7°, (8) 
and in the particle-hole conjugation, 

G,(r,r')^-7°G_e(r,r')7°, (9) 

both leaving the Dirac equation (jT]) invariant. Therefore, 
if boundary condition ([2]) deviates from the infinite mass 
confinement ([6|), our boundary problem exhibits no par- 
ity invariance and, in view of Eq. (I9|), no particle-hole 
symmetry. The symmetry breaking persists in the limit 
of the zigzag edge, Eq. ([5]). We therefore conclude that 
the parity breaking is due to the structural sublattice 
asymmetry of the zigzag- type lattice termination.^^ For 
practical calculations, we need to take into account de- 
viations of n± from z (e.g. due to potential-energy sub- 
lattice asymmetry^), because it eliminates the Green's 
function singularity at e = characteristic to dispersion- 
less zero-energy edge states. 

The connection between the parity breaking and the 
existence of the edge states can be established by explicit 
calculation of the Green's function from Eqs. ([T]) and 
The calculation details are given elsewherei^S Here we 
present the final result: 

o,..,. ^ (^)(^ + f^.) 

X (Gf,Uy,2/')/ + Gf;,(2/,y')S^) t^--ll^ (10) 



terms with the pole at e = hvrknxT- Let us examine, for 
instance, Eq. near the pole: 

GfJ,iy,y') - _Z!fi|(^a,e-i'=""i(«+^'). (13) 

s — nvrkuxT 

Clearly, the pole exists only if the unit step function 
Q{knzT) is nonzero, which determines the spectrum as 

e^-fe — hvTkrixT, kuzT > 0. (14) 

These equations are not yet restricted by the T symme- 
try. The T-symmetric spectrum follows from the con- 
dition that both equations (|14p are invariant under the 
simultaneous reversal of the chirality and wave-vector, 
T,k —T, —k. This imposes the following restrictions 
on iIt-: 

nxT = nx, n^r=Tnz, n = (n^;, 0, nj, = 1, (15) 

leaving a single free boundary parameter - the direction 
of the unit vector n. The edge-state spectrum is now 
manifestly Kramers degenerate and particle-hole asym- 
metric)^ 

e^f. = fivrkrix, rkn^ > 0. (16) 

The role of the parity breaking is quite apparent from 
the behavior of the edge-state Green's function (|13p : 
it vanishes identically for the infinite mass confinement 
(riz = 0) which preserves the parity symmetry (see, also 
Eq. jSl)). 

The knowledge of the spectrum (fT6|) is not sufficient 
to interpret the STS measurements p^i^^ as they provide 
information on the local DOS rather than dispersion Srk- 
In the next section we use the full Green's function ((TO)) to 
calculate the local DOS of the system. We will see that in 
addition to the exponentially localized states (fT3|) there 
is another type of edge states decaying algebraically as a 
consequence of the lack of the energy gap in the 2D bulk. 
This distinguishes our system from, e.g., topological in- 
sulators where bulk excitations are fully gappedi^i^i^^ 



III. LOCAL DOS AND PSEUDOSPIN 
POLARIZATION 



G 



(11 - ^ fp~qiy+y') _ p~q\v-y'\\ 



2(£ — hvTknxr) 



(3) , ,^ _ fc + gn^T - Teuxr/fiv 



(11) 



2(e — Uvrkuxr) 



-q{y+v') ^(12'^ 



where q = ^k"^ — /fflv"^ and k is the wave number. In 
Eqs. pT|) and (fT^ the edge states are described by the 



A. Particle-hole symmetry and the role of parity 

The spectral and pseudospin properties of the system 
are interrelated. Let us define the local DOS 

v±{e,v) - --^ImTr(/±75)Ge(r,r), (17) 

ZTT 

and local pseudospin polarizations 

P±(e, r) = - 2^IniTr (/ ± -i'')Y?G,{r, r) (18) 
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in terms of the RH and LH projections of the Green's 
function, G^. TunneHng spectra are determined by the 
total local DOS related to G(°) in Eq. (fTT|) : 



v{e, r) = + 



E ^^GfJ,{y,y). (19) 

T=±l,k 



Likewise, p+ +P- = — ;^ImTr J^^G^ is the net pseudospin 
polarization. It vanishes by T symmetry, since Eqs. (|15p 
yield p_ = — p+ . As a T-invariant characteristic of the 
pseudospin properties, we use the chiral pseudospin po- 
larization (GPP) related to G^^) in Eq. (fTS]) : 



p(e,r) =p+ — p_ = -ilmTr7^I]^Ge(r, r) = 

= -;^ E rImGf;,(y,y). (20) 

T=±l,k 

Integrating over k in Eqs. and (PO)) . we obtain 

2|e| 



^-^ ^ OvB "-I 



E 

T = ±l 



c?7 X 



(21) 



n^^ cos sin 7) + rn^rnzr sin 7 sin ( ^ sin 7) 







V / 


\\ 1 

\ '--^ 

w 


V A 




V ' 




11 m Zr 








y=6 nm 


E(eV) 






y=150 nm >\ 


-0.4-0.3-0.2-0.1 0.1 0.2 0.3 0.4 


-0.4-0.3-0.2-0.1 0.1 0.2 0.3 0.4 


C 1 




d 1 


Av 




«^ 8 =0.15 eV 


V 




0.25 


0.75- 










0.5- 










0.25 




-0.25 




10 


20 30 40 







FIG. 1: (Color online) Local density of states in units of 3/47r 
eV^^nm^^: (a) vs. energy for different at i/ = 4 nm; 
(b) vs. energy at different positions for — 0.95; (c) vs. 
position for opposite-sign energies and Uz — 0.95. Inset: local 
STS geometry and orientation of the unit vector n [Eq. (|15|l ] 
determining the boundary condition ([2}. (d) Asymmetric 
local DOS Av = //(e) — vi—e) vs. energy at t/ = 5 nm for 
riz — 0.95. The data are for ii = lO'' ms~^ and Ux < 0. 



p(e,,) = - Y: ^^7^^;^ dye -I "..I 



T = ±l 



E 

T = ±l 



(i7tan 7 x 



(22) 



rn2^ cos ( ^ sin 7) -I- n^;^ sin 7 sin ( 2gi sin 7) 
n^T- -I- tan^ 7 

It is now easy to see that the particle-hole symmetry 
is controlled by the parity. In the broken-parity state 
with n_|_ ^ n given by Eqs. (|15p . the summation over 
chiralities t = ±1 in Eq. (|2T|) yields an asymmetric DOS 
as a function of energy, e: 



2|e| 



29 (en^n^) 



'kTl'^v^ nv\n 



d'-f X 



(23) 



nl cos sin 7) +12^12^ sin 7 sin (2^ sin 7) 
X - 2 

+ tan 7 

The same is true for the GPP (l22l): 



2n2e(enj,n2) 
p e, y = — — ■ dye "-I --I 



2kl 



d-f ta.n 7 X 



(24) 



n2C0s(^sin7) -I- n^; sin 7 sm sm 7 j 
nl + tan^ 7 



For comparison, if the parity is preserved for n+ = n_, 
the DOS (|2T|) appears to be an even function of energy: 



«^(e,y) = 

_2\e\nl 

j^2f^2y2 



-dye- 



dj 



cos(^sm 



iin 7) 



tan^ 7 



(25) 



However, the requirements for the parity symmetry, 
n+ — n_, are incompatible with conditions psp for the 
T symmetry. The only exception is the infinite mass 
confinement limit Uz — > 0. In literaturCfii'^'^^i^l T- 
symmetry breaking on zigzag graphene edges has been 
discussed in connection with their possible intrinsic mag- 
netism. It is still unclear whether the T-symmetry break- 
ing in the boundary condition, Eq. ^ has anything to 
do with the edge magnetism. We will therefore limit our 
analysis to the T-symmetric case (|T5|) . 



B. Energy and position dependence of the local 
DOS: Analysis 

In Fig. [1] we plot the local DOS (|23|) as a function of 
energy e (in eV) and position y (in nm) for the Fermi 
velocity v — lO^ms^^. These units and parameters are 
typical for STS in graphene. Panel (a) shows an asym- 
metric peak, due to the edge states, emerging on top of 
the linear DOS as the boundary condition varies from 
the Berry-Mondragon type {riz = 0) to the zigzag type 
{\nz\ 1). In the latter case, the DOS ^ stiU fails to 
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FIG. 2: (Color online) Chiral pseudospin polarization (a) vs. 
position and (b) vs. energy for = 0.95 and < 0. 



recover the particle-hole symmetry because of the bro- 
ken parity [see, Eq. ([S])], which is formally described 
by the singular energy-dependent factor Q{enxnz). The 
crossover between the Berry-Mondragon and zigzag cases 
can in principle be induced by a staggered mean-field sub- 
lattice potential whose strength is parametrized by the 
angle between n and 

For 1 71,2 1 — > 1 the edge modes are much slower than 
the bulk ones, and have a small characteristic energy, 
w hv/2y X \nx/nz\. This can explain the observed DOS 
asymmetry on the scales of 20 — 50 mcV^^. For = 0.95, 
the peak position, e « 25 meV, and its overall behavior 
[panel (b)] agree very well with the observations (see, e.g.. 
Fig. 5 in Ref. 15). As we neglect possible level broaden- 
ing, the peak looks somewhat higher and narrower than 
in the experiment. Also, in agreement with the findings 
of Li et ali^, i/(e) approaches the symmetric Dirac DOS 
away from the edge. The position dependence of the 
DOS [panel (c)] shows that at the edge ly reaches either 
a maximum or a minimum depending on the presence 
or absence of the exponential term in Eq. (j23p . which is 
controlled only by the sign of e. The Dirac waves, inci- 
dent from the bulk and reflected from the edge, interfere 
yielding an oscillatory contribution (third term) in the 
DOS decaying as y^"^^^ on the scale of 10-20 nm. 

Panel (d) demonstrates similar oscillations in the energy 
dependence. 



IV. GLOBAL EDGE DOS AND PSEUDOSPIN 
POLARIZATION 

A. Relation between the edge DOS and pseudospin 
polarization 

Let us define the DOS and CPP of a finite region of 
space, < y < w, as the following dimensionless inte- 
grals. 



Ne{e,w) = hv dy [ iy{e,y) 



P{s,w) ^hv dyp{e,y). 
Jo 



, (26) 
(27) 



In the first equation we subtract the bulk Dirac DOS, 
so that Nf.{e,w) contains the contribution of the edge 
states only. By contrast to their local counterparts, it is 
convenient to call Ne{e,w) and P{e,w) the global edge 
DOS and global CPP, respectively. 

Inserting Eqs. ^ and ^ into Eqs. ^ and ^ 
and integrating over position y, we find 



"2 

X J d-y 



P{e) = 

"2 

X J dj 





26 (e n^n^) 



1 



(28) 



sin sin 7) 
sin 7 11^ + tan^ 7 



2nzQ (e n^riz) 



l^cos(^sin7) 
nl + tan^ 7 



2sgne 



(29) 



sm7j 



(%^sin7) 



sin 7 1 -|- nl cot 7 



1 + cot 7 



It is instructive to discuss first the limit w — » cx), when 
the integrals in Eqs. ((28| and ((29)) can be evaluated ana- 
lytically. In this case the integrals with the rapidly oscil- 
lating cosine function vanish, while those containing the 
sine function should be evaluated with care since the ratio 
sin (^i^ sin 7) /sin7 becomes singular, tt sgn(£) (5(sin7) 
as w 00. After integrating with the delta function 
(5(sin7), we have 



C. Local chiral pseudospin polarization 

To conclude the analysis of the local properties, in 
Fig. m we plot the CPP given by Eq. From Fig.[D(a) 
we see that the CPP has a purely boundary origin as it 
decays to zero in the bulk. Apart from the presence of 
the oscillations, both position and energy dependences 
of the local CPP differ significantly from the correspond- 
ing behaviors of the local DOS [cf. Figs. [2] and [T](c,d)]. 
Although not obvious in the local quantities iy{e, y) and 
p{e,y), in the next section we establish a direct relation 
between the appropriately defined global CPP and DOS. 



, , 2Q (e nxTiz) 2nj:nzSsne f d'y 
Ne{e) = — - 1 ^ ' ' 



+ tan^ 7 ' 



Pie) 



2nzQ [enxUz) 2na;Sgn(r 



c?7 



TT J 1 + nj, cot 7 




The remaining integrals are easy to evaluate^S.. The re- 
sults are 



1 -|- n^sgn(e nS) 
Ne{e) = , , — — - 1, 



(30) 
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FIG. 3: (Color online) Typical behavior of the global edge 
DOS, Ne and chiral pseudospin polarization, P as functions 
of the width, w (see, Eqs. ^ ~ ^ ). 



p^^^ ^ sgn(sn,) + ».^ ^ l + n sgn(sn.)^ ^^^^ 



From Eqs. ([50)1 and ([51]) we find the relation between 
and P, 



7Ve(e) = |P(e)| -1. 



(32) 



Being a scalar, does not depend on the sign of P 
which is reversed under the transformation n —^ — n. 

According to Eq. ([32]) . for the zigzag graphene edge 
rizl — » l{nx 0) the absolute value of the CPP becomes 
equal to the edge DOS: 



Neis) « \P{e)\ 



2Q{enxnz 



(33) 



This means that the CPP can, in principle, be detected 
through measurements of the global edge DOS. The lat- 
ter, in turn, can be probed by tunneling, as we discuss 
in the next subsection. Before going to that question, we 
wish to point out that the correlation between P and Ne 
exists for finite values of w as well. Figure [3] shows the 
functions Ne{w) [Eq. ^] and P{w) [Eq. clearly 
approaching the relation ([32]) for w > 100 nm. Note that 
the non-oscillatory components of P and obey the 
relation (1321) at much smaller w. 



B. Tunneling spectroscopy 

Our proposal for tunneling spectroscopy of the global 
edge DOS exploits the particle-hole asymmetric nonana- 
lytic energy dependence of iVe(e) [Eqs. ((30)) and ((33)) ]. It 
is essential that the particle-hole asymmetry persists in 
the case of zigzag-terminated graphene {\nz\ — > l,nx — > 
0) because this is an experimentally accessible system. 

It is known^^''^" that a strongly energy-dependent DOS 
reflects in the differential electric conductance, g{V) of 
a tunnel junction between the system of interest and a 
metal where the DOS is almost constant near the Fermi 




FIG. 4: (Color online) Suggested tunneling device for deter- 
mining the density of the edge states in graphene. A metallic 
film (M) is deposited on top of a zigzag-terminated graphene 
sheet (G) forming a strip-like lateral contact of width, w. The 
device resistance is assumed to be determined by the tunnel 
barrier (dark gray area) so that the voltage drop V predomi- 
nantly occurs between the overlapping parts of M and G. 



energy. Here we consider a lateral tunnel contact be- 
tween a zigzag-terminated graphene sheet and a metallic 
film, as shown in Fig. HI It is assumed that the voltage 
drop, V, occurs predominantly across the tunnel bar- 
rier in the contact area, which determines the junction 
resistance. Under such condition, the conductance can 
be calculated using the tunneling Hamiltonian approach, 
which is well described in the literature (e.g. Ref. 13(1 ), 
with the following result: 



l{V,T)^go deN{s) 



^df ie~eV,T) 
dieV) 



(34) 



Here N{e) is the DOS of graphene in the contact region 
(dark gray area in Fig. |4]), f{e — eV,T) is the Fermi- 
Dirac distribution of the tunneling quasiparticles at volt- 
age V and temperature T, and the constant go absorbs 
the energy-independent parameters of the metal and tun- 
nel barrier. As we are interested in the low energy regime 
y,T ^ 0, in Eq. ((34)) we can neglect inelastic tunneling 
processes (e.g. phonon emission) 

The DOS N{e) contains the contributions of both bulk 
and edge states. Since the bulk DOS is a symmetric 
function of energy (oc |£|), it can be eliminated by taking 
the difference: 



Ag(V^,r) = g(F,T)-g(-V^,T) 
= go 



de m-e) - N,{e)]^I^L_^^. (35) 



It contains only the particle-hole asymmetic edge DOS, 
Ne{s), given by Eq. ((28)) . where w coincides with the 
width of the lateral tunnel contact [see, Fig.U]. 

In the limit w oo, we use Eqs. ((31)) and ((32)) to 
evaluate the integral in Eq. ((35| : 



Ag(l/,r) = go Ap tanh 
Ap = |P|e>o - iJ'U^n = 2 



eV 



(36) 



71 

P\e<o = 2^. (37) 
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FIG. 5: (Color online) Zero-temperature conductance Ag = 
g{V) — g(— V') vs. voltage (in volts) for different widths of the 
tunnel contact, w [see, also Fig. |4]. 

The conductance asymmetry Ag(y, T) reflects the 
nonequiUbrium quasiparticle accumulation that builds 
up near the graphene edge in response to the current flow 
between the systems. For \eV\ > 2kBT the conductance 
Ag{V,T) saturates at ±goAp, where Ap is the differ- 
ence in the absolute values of the CPP for the positive- 
and negative-energy edge states. Such a nonlinear be- 
havior can be used to detect the edge state as well as the 
existence of the pseudospin polarization. At zero temper- 
ature r = the voltage dependence in Eq. ([36]) becomes 
singular (cx sgn(ey)). This is specific to the w = oo 
limit. As shown in Fig. [5l the singularity is smeared due 
to finiteness of the contact width, w, so that Ag(y, 0) 
saturates at voltages larger than the value oc w^^. The 
data in Fig. [5] are obtained by numerical integration of 
Eqs. dSSl) and (1251) at r ^ 0. 

V. SUMMARY AND DISCUSSION 

We have considered the boundary problem for 2D 
Dirac fermions, Eq. ([1]) and in which the time- 
reversal invariance is preserved at the expense of the 
Dirac fermion parity. Using the Green's function solu- 
tion, we have shown that the broken parity manifests it- 
self in the density of the edge states and their pseudospin 
polarization, both exhibiting a nonanalytic particle-hole 
asymmetry. The zigzag graphene edge with its inherent 
structural asymmetry is an example of the realization of 



the Dirac fermion parity breaking. Taking into account 
additionally the potential-energy sublattice asymmetry 
near the zigzag edge^, we obtain the local DOS consistent 
with the experimental data of Niimi et alJ^ We have also 
established a direct correspondence between the pseu- 
dospin polarization and the density of the edge states, 
and suggested how to detect them in a tunneling exper- 
iment. The proposal relies on the broken particle- hole 
symmetry resulting in an asymmetric nonlinear contri- 
bution to the conductance, Eq. (j36p . in a tunnel junction 
between zigzag-terminated graphene and a metallic film 
[see. Fig. g]. 

It is interesting to discuss possible implementations of 
the strong nonlinearity of the conductance, Eq. ([55)) . We 
suggest that it could be used for detecting weak electric 
signals and their polarity. The operation of such a device 
would exploit the two different states of the tunnel junc- 
tion corresponding to the conductance values at positive 
and negative bias voltages [see, also. Fig. [S]. Let us as- 
sume that the system is initially in one of these states. 
Then, under externally induced change in the bias volt- 
age the system can switch into the state with the other 
(lower or higher) value of the conductance. For the zigzag 
graphene edge (jn^l — > l,nx — > 0), the conductance dif- 
ference, Eqs. ([5S)) and ([57)) . is very significant and, there- 
fore, should be detectable. The progress in characteri- 
zation of graphene edges^^ may eventually lead to more 
understanding of how the interfaces needed to test our 
finding can be fabricated. 

We noticed that even though the conductance, Eq. ([551) 
vanishes for the Berry-Mondragon confinement^^ (jiz — > 
0), this case is still nontrivial, because the pseudospin 
polarization P{e) ([5^ is not zero: P{s) = sgn{enx)- 
This is confirmed by the numerical data for Uz = 0.05 
in Fig. [3] The nonvanishing P{e) comes from the oscil- 
latory (interference) term in Eq. ([^^ . decaying as y"^^"^ 
on distances 50-100 nm [see. Fig. j^fa)]. For a given chi- 
rality (valley), the long-range polarization implies viola- 
tion of the T symmetry on mesoscopic scales, which may 
have some connection to recent studies of the level statis- 
tics in graphene quantum dotsi^i^^ Also, such long-range 
polarization may coexist with magnetic correlations pre- 
dicted for zigzag graphene edges ^iSiiS^iSl since they are 
expected to decay on much shorter distances<21 
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